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Abstract 



The spontaneous magnetization in Chern-Simons QED% is dis- 
cussed in a finite temperature system. The thermodynamical poten- 
tial is analyzed within the weak field approximation and in the fermion 
massless limit. We find that there is a linear term with respect to the 
magnetic field with a negative coefficient at any finite temperature. 
This implies that the spontaneous magnetic field does not vanish even 
at high temperature. In addition, we examine the photon spectrum 
in the system. We find that the bare Chern-Simons coefficient is can- 
celled by the radiative effects. The photons then become topologically 
massless according to the magnetization, though they are massive by 
finite temperature effects. Thus the magnetic field is a long-range 
force without the screening even at high temperature. 
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1 Introduction 



The 2+1 dimensional gauge field theories with Chern-Simons (C-S) term 
have had a lot of motivation to study. In addition to their own theoretical 
interests, they are expected to describe the physics of the planar systems like 
the quantum Hall effect |l|, ^ and the high T c superconductivity || . Many 
interesting and remarkable properties of these theories have been found in 
past a decade, such as their chiral dynamics [§], ||, the spontaneous parity 
violation by radiatively induced C-S terms 0. A few years ago, it was 
also found that in a model of three dimensional quantum electrodynamics 
(QED 3 ) with a bare C-S term the spontaneous magnetic field can be stable 
by the fermion loop effects, which eventually triggers the dynamical Lorentz 
symmetry breaking (Hosotani pj). He has calculated the energy densities as 
a function of the magnetic field B in the fermion one-loop level and showed 
that the minimum of the energy density can be at B ^ under a consistency 
condition for non-zero magnetic fields. It has been found that this condition is 
also regarded as the sufficient condition for the spontaneous generation of the 
magnetic field as well ||. According to the spontaneous Lorentz symmetry 
breaking, the physical photon, which is topologically massive in the tree level 
due to the bare C-S term, becomes massless as a Nambu-Goldstone boson. 
This mechanism of the spontaneous magnetization has been extended into 
the cases with massive fermions with non-zero fermion density M. Even 
in the case of no bare C-S term, the spontaneous magnetization can occur 
by adding heavy fermions because the C-S term induced by them can play 



the same role as the bare C-S term [10]. The spontaneous magnetic field is 



expected to be important in the connection with the problem of the fractional 
quantum Hall effect 0]. 

In this paper, the thermodynamic properties of the spontaneous mag- 
netization in QED 3 with a bare C-S term are investigated. Some of the 
thermodynamic properties of the model without external fields have been 
studied to some extent in the different contexts so far. The specific heat for 
different values of the statistical angle is analyzed in ref. [fLlJj . The restora- 



tion of the parity at high temperature limit has been discussed in ref. |12| . 
in which it has been found that the radiatively induced parity violating C-S 
term vanishes only in the high temperature limit. Our aim here is to exam- 
ine whether the spontaneous magnetic field, which is dynamically generated 
at T = 0, vanishes or not in a high temperature region. The model we are 
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studying here is based on ref. M, in which all fermions are very light so 
that a non-zero bare C-S term is needed in the Lagrangian for the magneti- 



zation [10]. In the finite temperature system, the energy density is replaced 



by the thermodynamical potential. The quantum fluctuation of the thermo- 
dynamical potential is investigated in the fermion one-loop level here. The 
coefficient of the linear term with respect to the magnetic field in the poten- 
tial is calculated as a function of temperature as well as the classical part 
of the potential. Negative values of the coefficient mean the existence of the 
spontaneous magnetic field at finite temperature. We find that the coefficient 
becomes more (negative) large behaving like ~ — log T in high temperature 
region. Thus the spontaneous magnetic field does not vanish even at high 
temperature. 

In addition, the spectrum of photons is also investigated in the finite tem- 
perature system. At zero temperature, the photons are topologically massive 
at tree level. The mass amounts to the C-S coefficient. If the magnetiza- 
tion occurs, the bare coefficient of the C-S term is exactly cancelled by the 



radiative effects |L4], so that the photons becomes massless. Therefore 
the photons are regarded as the Nambu-Goldstone bosons of the dynamical 
breakdown of the Lorentz symmetry |7|, ||. On the other hand, at finite tem- 
perature, the Lorentz symmetry is explicitly broken by thermal effects. The 
photons are then no longer massless even if the magnetization occurs. In 



general, they eventually have the thermal masses (T6|, [L7|, |18j. The magnetic 
part of the field strength, however, become statically massless and the mag- 
netic force is a long-range force without the screening [[TJ] if the topological 
mass is cancelled out. We find that, in our model with massless fermions, the 
bare topological mass is completely cancelled by the radiative effects even at 
finite temperature if the magnetization occurs. This implies that the static 
mass of magnetic fields vanishes according to the magnetization even at high 
temperature. This is consistent with the result of the non-vanishing mag- 
netic field by the thermodynamic potential approach. It is also found that 
in the case of massive fermions the topological mass is cancelled out only in 
the 'chirally symmetric' cases with massive fermions. 

Furthermore we naively realize the reason of the spontaneous magnetiza- 



3 In our previous preprint |T^| , some insufficient treatment has been done especially in 
the evaluation of the contribution of the zero-mode to the coefficient. These are correctively 
changed in the first half of this paper. 
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tion from the viewpoint of the difference of the free energies of free photons 
between B and B = 0, where B is the magnetic field. It is found from 
rough estimation that the linear \B\ term with a negative coefficient appears 
in the difference. This explains the magnetization at finite temperature to 
some extent. These characteristics of the magnetization at finite temperature 
might be interesting in the connection with the physics such as the high Tq 
superconductivity. 

This paper is organized as follows. In the next section, some general 
properties and statistic characteristics of the model will be presented. In 
Sec 3, the spontaneous magnetization in this model in T = case will be 
reviewed briefly for a preparation for T ^ cases in order to clarify the 
essential point for the magnetization and also to fix some notations. The 
strategy to extend the investigation to the finite temperature system will be 
explained then. Sec. 4 will be devoted to the calculation of self-energies of 
gauge bosons in the one-loop level and the weak-field approximation. In Sec. 
5, the results that the spontaneous magnetic field does not vanish will be 
presented in the fermion massless limit by analyzing the effective potential 
both analytically and numerically. In Sec. 6, the spectrum of the physical 
photons will be examined in the finite temperature system. In Sec. 7, the 
reason why the spontaneous magnetic field does not vanish will be discussed. 
Also some comments on our results will be in order. In the last section, 
we will summarize our results. Some definitions, details of derivations and 
formulas are summarized in Appendices. 



2 Model in Finite Temperature System 

We consider the model of QED 3 with a bare Chern-Simons term de- 
scribed by the Lagrangian 

C = -~F^F^ - ^A,d v A p + J2$ a m*d» + q a A,) - m a } ^ a , (1) 

where gamma matrices are defined as 7^ = (f] a o' 3 , icr 1 , iu 2 ). There can be two 
types of photons according to the choice of r\ a = (i/2)tr7°7^7^ = ±1. Since 
the model ([!]) has charge conjugation invariance, we can take the electric 
charges q a to be positive without loss of generality. Also since the trans- 
formation m a <-> — m a is equivalent to r] a <-> — rj a , we can consider m a as 
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non- negative. We then call the fermion with 7] a = ±1 as ^-fermion, respec- 
tively. 

Under the existence of a classical magnetic field B ((0|v4 M |0) = —5 IMl x 2 B), 
the energy spectrum for the rj + fermion takes the form Jl9| 



E* = e(q a B) x u a , 
E a = ± w« (n > 1) 



(2) 



where e(x) is the usual sign function, lu® = ym 2 + 2n/Z 2 and l\ = l/\q a B\. 
Then the general solution of the Dirac equation is 

where u c np = U c u*, U c = 72 and the same for w^ p . The positive energy 
solution u np and the negative energy solution w np in eq. (§) are 

„„„(*) = J= e - W) f "»«) j , (4) 



w np [x) = —f=e 



where we dropped the fermion index a for simplicity, and k = 2np/Li (p: 
integer), f = x 2 // - Z/c and u n (f) = (-l) n 7r- 1 / 4 (2 n n!)- 1 / 2 e« 2 / 2 rf n /^ n e" e2 - 
The solution for r?_ fermion is obtained by making the transformation (t — > 
— £) from corresponding ?7 + cases. There is an asymmetry between positive 
and negative energy solutions in the lowest Landau levels. The vacuum 
expectation value of the total charge is (see eq. (|A.2|j ) 

<0|Q|0> = E?.*(B)(".4)> ( ? ) 

where the number of degeneracy is given by = J2 P 1 — \laB\ /2n x / d 2 x 
and we introduced the filling factor v a = (0\u a \0) for each fermion ip a , 

= f E P <a 0p /iV s a , (rj a e(B) > 0), 
a 1 E P &U*W. ( Va e(B)<0). 1 ; 
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We here consider the completely filled (u a = 1) or the empty [y a = 0) cases 
for zero temperature. From the time component of the equation of motion 
with the charge, we have a relation / d 2 xnB = (0|Q|0). The consistency 
condition for B ^ is then obtained from this relation and eq. (|7|) as 

In this paper, we would like to study this model in the finite temperature 
system. The partition function for the fermion sector is given by 

Z = tre-^-^ 

a 

X f[ (l + e -0( w S-A 4 »9»)) iV? + e -i3«+»aq a )^ N ° e /3<iV s a ; ( 1Q ) 
n=l 



where (3 = 1/T and nQ = Y.al^aQa (see eqs. ( |A.1|) and ( |A.2|) ). The total 



charge Q = J2 a Qa in the finite temperature system is obtained by using (|I^ 



Note that our definition of the chemical potential is a little different from 
the ordinary one and q a ^ a in our notation is normally called the chemical 
potential in the literature. The condition for B ^ in the finite temperature 
system is given from the time component of the equation of motion as k = 
Q/(VB), where V is the two-dimensional volume. In the massless limit 
{m a — > 0), this becomes by using eqs. (|T0l) and (|TT|) 

k = e(B)£||{±taiihQ^) 

+ ^ ( e /3«-9a/^) ZjZ 1 _ e /3«+9aMa) + \) ' 

The chemical potential fi a is determined as a function of T and B by eq. 
(|T2| ) or by the following relations which are induced from eqs. (||) and (12): 



tanh (-q a (3\n a \ 

00 / ] i \ 

+2 =i Ve^<-^l^l) + 1 ~ e /3«+g a |Ma|) _|_ ij ' ( 13 ) 
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Figure 1: The chemical potential \fi a \ in the fermion massless limit. By 
assuming the consistency condition (12), fi become a function of B and T. 
It can be seen that \y, a \ is a monotonous decreasing function of T with the 
definite magnetic field \B\. The solid line is the case of \q a B\/ h 2 = 1/2. The 
dashed line is the case of \q a B\/n 2 = 1/4. 

and e(/x a ) = 2e(B)rj a (i , a — 1/2). The asymptotic behavior of \fi a \ is obtained 
from (|i~3|). We can easily see that \fi a \ ~ cj"/2 = ^\q a B\/{2q 2 ) in the low 
temperature limit (/? — > oo). At high temperature, we also see that \fi a \ 
behaves like ~ /3\B\. We show the behavior of the chemical potential as a 
function of T at some values of B (Figure 1). It can be seen in the figure 
that \fi a \ is a monotonously decreasing (increasing) function of T for fixed 
B > (B < 0) and is also a increasing function of \B\ with fixed T. We also 
find from the relation (O) that for small /3 2 \B\, 

l/ia| ^41u2 /9|jB| + ° (|jB|2) - (14) 

Note that the non-zero fi a here comes from the asymmetry of the lowest 
Landau levels. In the case of B = 0, the asymmetry vanishes so that fi a 
becomes zero. 
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Figure 2: The thermal behavior of the filling factors (P ) for fixed \B\. The 
solid line shows the behavior of (y + ). The dashed line is that of We 
here set v + = 1 and v~ = at T = 0, taking account of eq. (9). 

Now we comment on the thermal behavior of the filling factors. The choice 
of the filling factors at zero temperature is taken to satisfy the consistency 
condition (|9|) with k / 0. The thermal average of filling factors is easily 
calculated by using (§) as 



When we set v + and v~ (the filling factors for rj + and ?7_-fermions, respec- 
tively) into 1 and at T = respectively, the thermal average (£> + ) ((£")) be- 
comes the monotonous decreasing (increasing) function of T and (u^) —>■ 1/2 
at large T (see Figure 2). This merely shows that the probability that the 
?7-|--fermion exits in a lowest Landau level is multiplied by no statistical weight 
at high temperature limit. 
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(15) 



e -/3q a [i a ri a e(B) _|_ ]_ ' 
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3 Spontaneous Magnetization, Strategy to Fi- 
nite Temperature System 

In this section, we will show our strategy to study the spontaneous 
magnetization in the finite temperature system. To this aim, we start from 
a short review for the zero temperature case, which have been discussed at 
first by Hosotani Next we shall discuss the extension to the finite 

temperature case. 

At zero temperature, Hosotani has shown in the same model that the 
spontaneous magnetic field can be stable and then triggers the spontaneous 
breakdown of the Lorentz invariance. The effective potential has been calcu- 
lated in a weak field approximation and the fermion massless limit J7|. The 
classical part (A£^ in ref. @) of the potential yields the contributions of 
the orders of |I?| 3//2 (matter) and \B\ 2 (photon) with positive coefficients. 
The minimum of the potential of the classical part is then at B = 0. The 
radiative correction to the potential drastically changes the situation. The 
deviation of the fermion-loop contribution to the effective potential between 
B and B = 0, £ rad (= A£ (1) - AS^, in ref. @), is given by 

.rad_ _£ f d 3 P ^ (l+Il Q ){l + (j^l Q -fU 2 )/p i }-(K-U 1 )yp i 

L ~ 2 J (2vr)3 111 (fl - 0) ' 1 J 

where n i (= Sail"), (i = 0, 1,2) are defined from the one-loop contribution 
of the vacuum polarization tensor: 

r^(p) = (p^-p 2 g^)U {p)+ie^p p U 1 (p) 

+(i _ _ 8 M ){^p v -fsn^ip) - n„(p)). (17) 

In the case of Vg a = e(> 0) and = Nj = N (called the chirally symmetric 
case), putting the filling factors v + = 1 and v~ = on taking account of 
the condition @, £ rad is then calculated in the fermion one-loop level by 
taking a weak field approximation and the fermion massless limit: £ rad = 

—en/if 2 tan _1 (4/7r) \B\-\ . This implies that there cannot be the minimum 

value of the effective potential at the point of B = so that the spontaneous 
magnetization occurs. Then by the relation (0\[iM 0t , F°i (x)]\0) = e^B ^ 
where M^ v are the generators of the Lorentz transformation, it is found that 
the non-vanishing magnetic field B induces the spontaneous breakdown of 
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the Lorentz invariance. Note that the essential point for this result is the 
existence of the non- vanishing negative coefficient of the linear \B\ term in 
the effective potential, which comes from the radiative correction part £ rad . 
If we set the bare C-S parameter k into zero, the linear term vanishes and 
the symmetry is restored. 

From now, we discuss the extension of the effective potential to the finite 
temperature system. The effective potential (energy density), £(B), is then 
replaced by the thermodynamical potential Q(B,T,fi) POfl . In the case of 
B = 0, Bralic et. al. has calculated the potential in this model in the 
different context fll |. In the B ^ case, the classical part of the fermion 



sector is given by 

Q^ = -±\nZ, (18) 

where V is the two-dimensional volume and Z is defined in eq. (|H]). We 
show an expression of Q f Q rmi in Appendix 1. We confirm that there is no 
linear B term in fig" there p6| |. The quantum fluctuation part of Q is 



more important here because we are interested in the spontaneous magnetic 
field for which the quantum part is essential in the case of T = 0. The 
extension of this part to the finite temperature system is well performed by 
the Matsubara method The fluctuation interaction part, £T ad , of the 

thermodynamical potential is then obtained from the eq. ( |16D just by taking 
the replacement 

p — > ip 3 = i — — , (m : integer), (19) 

(IP % r ^ . . 

£ - (20) 

and Tl°;(p 2 ,B) — > II^p 2 , B, (3) m , (i = 0,1,2) for bosonic loop integrations. 
As to the fermion loops included in IT, we employ the replacement with 
imposing anti-periodic boundary condition, 

k -> ik 3 = % h q a ^ a , {n : integer), (21) 

— ^ ^E- (22) 

27T £V 
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Thus the deviation of the quantum fluctuation part of the potential is for- 
mally expressed in the finite temperature system by 



n 



rail 



2(3 



E 



d 2 p 



x In 



(i + flo) {1 + (^flo + ?n 2 )/(rt + p 2 )} + (« - noVbi + p ; 



(5-0) 



(23) 



We, however, will not treat this expression directly. Since the non-zero 
coefficient of the linear B term is essential for the spontaneous magnetization, 
what we have only to know in detail for our purpose is just the coefficient 
C(/3,/i ) m^ rad , * 



n rad = C(P,no)\B\ + (higher orders of \B\), 



(24) 



where /x = h(/3)\b=o is zero in the fermion massless limit because of eq. 
flUf) . We will concentrate our attention into C{0) = C(f3,0) from now. The 
coefficient is decomposed into two parts as 



C(/3) 



rad 



OB 



+ 



B=0,n=0 



d/i a dB 



B=0,At=0 



(25) 



The first term is obtained from eq. ( p3|) as 



CM 



- E 



d 2 p 



dU 
dB 



{pl+p 2 + {plti +p 2 ti 2 )} 



+(i + n ) \pI 



(I ^2 5II 2 
P 



dB * 3B 
1 + Uo) {(£ + ?)+ 
p 2 3 U + fU 2 } + (K-U 



2(« — IIi 



glli 
' ~dB 



B=0,n=0 



(26) 



The second term in eq. fl2"5] ) is the contribution of the chemical potential 
fi. We however easily see that C^{(3) does not contribute. Since n rad is the 
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deviation of the quantum fluctuation between at B ^ and at B = 0, we 
can easily see that C^ft) becomes zero by virtue of eq. fll3|). 

As long as the coefficient C(/3) is negative, the spontaneous magnetic field 
does not vanish. In order to evaluate the behavior of C(0) by using (p6|) , we 
need to calculate IT's. In the next section, we will show the details of the 
calculation of HVs. Since we are interested in the coefficient of the linear \B\ 
term, Cq(/3), the calculation will be performed up to the order 0(B l ) there. 



4 Boson Self-energies at finite temperature 

This section will be devoted to the calculation of the boson self-energies 
Ui(i = 0,1,2) at T ^ 0, which are composed of the fermion one-loop di- 
agrams. In the T = field theory, the fermion propagator S(x,y) in the 
classical field B is expressed by the proper-time method 0. We just quote 
the expressions and some character of S(x,y) in Appendix 2. The fermion 
one-loop contributions to the vacuum polarization tensors are then calculated 
by 

—tr{rSS(k)YS a (k-p)}, (27) 

where propagator Sq(p) is defined in eq. (|A.7|) and T Ml/ (= J2a^ afiu ) are 
related to £ a u{) in eq. (|T7j). 

Since we are interested in the coefficient of \B\ in f2 rad , we here take 
the weak field approximation up to the order 0(B l ) in the calculation with 
setting fi a = 0. (In Sec. 6, we will adopt a little better approximation to 
evaluate the spectrum of the photons.) We here show it only in the case of 
r7 + l-fermion in q a B > because other cases are connected with this case by 
eqs. ( |A.6|) and (|A.9|) . At first, we consider the contributions of fermions with 



v a = 0. We drop the fermion index a for a while just for avoiding complexity. 
The expansion of SqIp^Zq 1 by l// 2 = \qB\ becomes 



pZ _ m z _|_ le _ m z _|_ %t y 

= 5 ( P )^ 1(0) + 5 (p)^ 1(1) r 2 + o(r 4 ). (2* 



We then obtain the expansion of T^ u from eqs. (|27| ) and (|28| ) in the case of 
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0: V^{p)lZt l = r^(p)lZt m + T^(p) v u Z+ m l~ 2 + 0(/~ 4 ), where 



r""(p)£ 



% q 



% q 



d 3 k 
(2^ 

d 3 k 
(2^)3 



tr 



tr 



rS (ktloYSo(k-ptlo}i (29) 
^So(k)V ol »So(k-p) { °io]}- (30) 



+tr 



Since we can easily check T >lu (p) 7 l = Q 1 ^ to be zero due to the properties 
of the Dirac 7-matrices, we can regard (p^Z^ 1 as V^" '(p)^=J in this 
approximation. Next, let us consider the case of v = 1. The propaga- 
tor in this case is given in ( |A.11|) . It is convenient to see the deviation 
ST^ip)^ 1 = T^ipYZt 1 - r^0)2=o\ where 



ST^(py= +1 



1 q 



+i q 



d 3 k 



{tr [^S {kfZtl v f{k-p)] 
+tr[Yf(k)rSo(k-p)lZ+ 1 }} 



d 3 k 



tr Wf(kWf(k-p)}. 



(31) 



Since f(p) includes the part e 



-p 2 i 2 



(see Appendix 2), it cannot be taken the 



Taylor expansion by l/l 2 in a usual sense. We have to keep this part in our 
calculation whereas we drop terms with the part l~ 2 x e~ p 1 . The last term of 



r.h.s. in eq. (0) is then dropped in this sense. Hence, in our approximation, 
Ilj's are calculated by the following expressions: 



1 q 



d 3 k 



rtr 



(2tt) 3 

sr^ipyzt 1 = T^{p)i=J l + iq 2 

+ tr 



YS (k)^ YSo(k - p)S°2 ] + 0(l~ 4 ), (32) 
d 3 k 



(2tt) 3 

^f{k)YS Q {k~ P )T = t m 



{tr[rS (ky u Zt m Yf(k-p) 

2-l 2 p 2 \ 



o(r z e 



(33) 



Before proceeding to the calculation of IL;, we mention the consistency of 
our approximation scheme. The self-energies at T = (Ilj's) are calculated 
by using the expressions in eqs. (W2T) and (RBf) up to the order Oi^B 1 ). They 
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are presented in eqs. ( |A.12| ) ~ ( |A.17p in Appendix 3. We can see that 
these expressions are consistent with the calculation by Hosotani [[?[. His 
calculation of the self-energies by virtue of the proper-time representation 
includes the ultra-violet divergence which has to be renormalized. On the 
other hand, we do not have any ultra-violet divergence in our calculation. 
This is because the divergence part appears only in the term higher order 
than (9(|I?|), which we are now neglecting. 

Now we go to the case in the finite temperature system. By applying 
the replacements (|2l"l) and (^) (and p — > ip 3 ) to HVs calculated in eqs. 
(TATT21) ~ flA~TFD with eq. QOg) , we obtain n , fli and fl 2 in the weak field 
approximation, which are presented in Appendix 4. From now we take the 
fermion massless limit (m a — > 0). We also restore the fermion index a in the 
expression. For the case of v a = 0, Ilj's become then 



fl a , u=0 (p3,p\B,f3) 



i6 (p + ply/ 2 

+ O(\B\ 2 ), 



+ A a (p 3 ,\p\,P) + O(\B\ 2 ), (34) 



(35) 



2 

■^fl^ u=Q (p 3 ,p 2 ,B,P) 



Qa 



I6p 



if+pD 1 ' 2 - ^(pa, \p\,p) + o(\B\ 2 ),m 



p- 



where Aq and A\ are contributions of finite temperature effects from the 
fermion loops and take the forms as 



A a (p 3 ,\p\,P) 



2nfr Jo 



dk- 



fik + i 



(f+P 2 - 4k 2 ) 2 + lQk 2 p 2 3 + p 2 +pj~ Ak 2 



A a 2 ( P3 ,\p\,P) 



phi 



2np 2 

ql (p 2 + pi) 1/2 



dk 



2{p 2 +p 2 ) 
1 



i \ 

2 



(37) 



47T 



p- 



e /3k + I 
1 



dk 



e f3k + 1 



{p 2 +pj-4k 2 ) 2 + \Qk 2 pl 



(38) 
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where 



(*) = V2\(4k 2 -pl) 



{p 2 +Pl- 4k 2 ) 2 + 16k 2 p 2 3 + p 2 +p\- 4k 2 

1/2 



nl/2 



(P 2 +Pl~ 4A; 2 ) 2 + 16Jfe 2 p§ - (p 2 +p\- Ak 2 



We can see from eqs. fl37|) and (^) that I1q and IlJj; grow like ~ T at high 
temperature for fixed p%. This behavior has already been pointed out in ref. 



12 , where the calculation has been based on the real-time formalism. As to 



the case v — 1, we obtain 

fi a , u= M,p 2 ,B,p) 

fl a hu=1 (P3,P 2 ,B,(3) 



fl a , u=0 (p 3 ,p 2 ,B } f3) } 
Va q -±^-^\B\ + 0(e-v 2 /^) 

7T p z + £>3 



^lu=M,P 2 ,B,f3) = fll u=Q (p 3 ,p 2 ,B,fJ). 



(39) 
(40) 
(41) 



Note that only Ul u=1 depends on the magnetic field in this approximation. 
We can neglect the term of the order 0(e 



-p 2 /(Qa\B\)^ 



m eq. 



for the 



evaluation of C((3) because this term is sufficiently smaller than the linear 
term for small \B\. 

Therefore, we have calculated all the self-energies II j up to the order 
0{B 1 ) in the fermion massless limit. In the next section, we will investigate 
the character of the effective potential f2 rad by analyzing the coefficient of 
the linear term C{f3) by using the results in ( |3"3| ) ~ (|4"T|). 



5 Thermo dynamical Potential 

We here examine the effective potential in the finite temperature system. 
The purpose is, as we mentioned before, to decide whether the spontaneous 
magnetic field, that takes non-zero value at T = 0, vanishes at high tempera- 
ture or not. We shall analyze the coefficient of the linear \B\ term C$(f3) in the 
potential Q(B, /3, /i a ) by using Ilj's calculated in the massless limit in the last 
section. In the evaluation, we set q a = e > 0, and assume N~}~ = Nj = N. 
The filling factors v + and V- are set their value into 1 and at T = 0, respec- 
tively. The condition (|9|) becomes then k = e 2 N/(2ii). By using eqs. (|34]) ~ 
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as well as this condition, the coefficient Cq(/3) in eq. ( p6|) is expressed 



as 



W) = -^ m EJ o ^^y, (42) 

where a = 2nT / k and x = 

The function F m (x , a) in r.h.s. of eq. (|42|) is defined as 



777. ( X , ) 



x 2 (x 2 + (am) 2 ) + —x 2 \jx 2 + (am) 2 + v4 (x, a, m) 



4 

x x 2 [x 2 + (am) 2 j + —x 2 ^Jx 2 + (am) 2 — A 2 (x, a, m) 
+x 4 (x 2 + (am) 2 ) , (43) 

where 

A (x,a,m) = ^2 {^am, kx, 27r/(/ta)) x 2 {x 2 + (am) 2 ) , (44) 



A 2 (x, a, m) = ^2 (nam, kx, 27r/(/ta)) ( — ) . (45) 



x N 2 



Before the numerical analyses, let us see the asymptotic behavior of Cq((3) 
described in ([42[ ). At first, we consider the low temperature limit. Since A 
and A 2 have the factor l/(e l3k + 1) (see eqs. (|37|) and (|38D), these terms are 
rapidly dumped. Then F m (x , a) is replaced by x 4 x /(x , am), where 

f(x , am) = (x 2 + (am) 2 ) jl + + ^/x 2 + («m) 2 J j . (46) 

The coefficient is then calculated by making use of the Euler-Mclaurin's 
mathematical formula (see Appendix 6.) as below, 

PK I °° /*00 rv> f>OQ rp I 

Co(/3) - -^ 2^/ ds - , - / dx- ' 



2tt 2 [ ^ 7 f(x,am) Jo f(x,0)j 

eK [2 /» , f» , x ,. /- 00 , x _ . . 

a < — / ay / ax—- + lim / ax—- + O(a) 



2n 2 [a Jo Jo f(x,y) m ^°°Jo f(x, am) 
en f°° , 1 



7T 2 Jo 



poo I 

/ dr 2 + (9(a 2 ), (x = r COS6 1 , ?/ = r sin 6 1 ) 

Jo 1 + (f + r) 

fean-^ + OCT 2 )), (T~0). (47) 
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Figure 3: The behavior of the coefficient C{0) (/3 = 1/T) is shown by the 
solid line. It can be seen that the coefficient is negative and does not vanish 
even at high temperature. Rather, the absolute value becomes larger as 
~ logT. The point A which is the low temperature limit takes the value 
ne/n 2 tan _1 (4/7r). The dotted line is the contribution of the zero-mode. The 
dashed line is the summation of the other modes. 

Thus we indeed reproduce the result in [Tj. Secondly, at high temperature 
region, the contribution of the zero mode (m = 0) becomes dominant because 
that of the non-zero modes (m ^ 0) is relatively suppressed by 1/T as 



The contribution of the zero mode does not vanish in the limit and behaves 




24ttT' 



(T - oo). 



(48) 
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like 



Co(/3)m=0mode -> - 111 — , (T -f Oo). (49) 

Aixml k 

Eqs. (|48D and fl49|) show that Cq{0) does not vanish at any finite temperature. 
The numerical analyses show that the coefficient Cq(/3) is always negative 
at any temperature (Figure 3). The non- vanishing Cq(/3) means that the 
thermodynamical potential Q does not have its minimum value at B = at 
any temperature. 



6 Photon Spectrum at Finite T and \i 

At T = and /i = case the spontaneous magnetic field implies the 
dynamical Lorentz symmetry breaking. The massless pole then appears in 
the spectrum of the physical photons by virtue of the Nambu-Goldstone 
theorem. Namely, there is no spontaneous magnetization at tree level and 
photons then have a topological mass \k\ whereas the magnetization occurs at 
one loop level and photons then become massless. The massless-ness comes 
from the exact cancellation of the bare topological mass k by IIi(O) 0, ^ [To - ]. 
It has been well known that the cancellation is exact for all orders of the 
perturbation |L4]]. 

In the finite temperature system, the Lorentz symmetry is explicitly bro- 
ken by the temperature effects. On the other hand, we have just shown by 
the thermodynamical potential approach that the spontaneous magnetic field 
does not vanish at any finite temperature. It is of very interest to examine 
the photon spectrum in the finite temperature system. The photons then 
becomes massive by temperature effects even if the topological mass is can- 
celled out at loop levels. From the Lagrangian ([!]) the photon spectrum in 
the finite temperature system is determined by [112], pl| 



(p 2 -fl 2 L )(p 2 -fl 2 T )-p 2 fl 2 opo = 0, (50) 

where IT^ and Ht becomes the thermal masses for the longitudinal and trans- 
verse modes respectively when the topological mass n topo is zero. In our 
model U.l,t an d n topo are 

MP 2 ) = -P 2 Y.^o(P% (51) 
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Mp 2 ) = -^EnS(p 2 )+p 2 E^(P 2 ), (52) 

a a 

n toP o(p 2 ) = K-Err^ 2 ), (53) 



where we wrote p% = —ipo- By using eqs. (|34|) , (pq), fl39|) and (|4T|) , the 
expressions of IT^ and ITt at high temperature are obtained in the chirally 
symmetric 

fli = Uf ^{j- 1 )- (54) 

fl T = 4JV^P|( po _p), (55) 

where u 2 = (e 2 ln2/47r)T is the plasma frequency. These are consistent with 
the results by Klein-Kreisler et. al. [[L2| . 



At first, let us consider the case of n topo (0) = 0. In terms of the field 
strength, and IT-p are the thermal masses of electric fields and magnetic 
fields respectively. In the static limit (p^ = 0, p — > 0) ||17|| , we easily obtain 
= AN u p and = 0. The electric fields are then massive and become 
short-range force and thus screened, whereas the magnetic fields are massless 
and remains long-range forces. In dynamic cases such as plasma oscillation, 
both fields are massive with the same mass ~ u p . Secondly, in the case 
of ritopo(O) = k 7^ 0, the electric and magnetic fields become both massive. 



They have the same static masses ~ \jANUp + n 2 . Also, the dynamic masses 



of the photons become approximately ~ ^j2Nujp + k 2 . Therefore if R is zero 
at any temperature, we can obtain a certain relationship at finite temperature 
between the magnetization and the statement that the magnetic fields remain 
long-range forces. 

In the following we shall show that in our model, the topological mass R 
vanishes at any temperature: 

« = «-£)n;(o) = o. (56) 



We stress that this is not contradictory with the results by Lykken et. al. |24 
In our model, we at first put the fermion distribution at zero temperature. 
The chemical potential is then not a free parameter but a function of B and 
T, which is determined by the consistency condition at finite temperature. 
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Now we calculate Ili(O) in the weak field approximation, we here leave the 
fermion massive (but not so heavy). For the u = fermions, the contribution 
is easily calculated from ( |A.25|) by putting P3,p = and using the formula 
<HKM ) as 

fif °(0) = f— " ( 1 - a (m „_ n ,.\ I -, - „flf m „4-„„,,.„1 I 1 ) ' ( 57 ) 



47]- Y e P{m a -q a Ha) _|_ X e P{m a +q a n a ) _|_ J 

We can easily see that, when m a = 0, I1^ =0 (0) becomes zero at any tem- 
perature and density as in eq. (|35|). In the low temperature limit, eq. (|57|) 
reproduces the result in || 0. At high temperature limit, it becomes zero. 
Next, we consider the case of v = 1 fermions. We cannot use eq. (^) to take 
the limit of P3,p = because the approximation for eq. (|40|) is of the order 
(9(|i?|) which is not appropriate here. So we go back to the expression in 
eq. (|A.16|) . Here we consider IT" i/=1 (0) as the form getting from applying the 
replacement (|T^) , (^T|) and ( P2"D to eq. ( |A.16| ) . By using the formula ( |A.34| ) , 
we can decompose 5Ilf(= Ul v=1 — n^ =0 ) as 

<jft;(p) = <m?(p) + n?(p, (3) Theimo , (58) 

where 511" is the deviation between v = and z/ = 1 of the induced C-S term 
at T = 0, which has been calculated by using the proper time representation 
as 0, 

AirM = i - i - i (m) 

In eq. (0), when we neglect the contribution of higher orders of \B\ than 
the first order, the result in eq. ( |A.22| ) is of course reproduced. This becomes 
r] a ql/(27i) for Pfl = 0. The second term in r.h.s. in eq. ( |58"D is defined in eq 
( |A.30| ). We obtain the expression as (see Appendix 5.) 

> fJ'a) Thermo 

= VaqXpoiPo + 2m a )e- po ^ +2m ^ {6(p ) - 9(- Po - 2m a )} 
j 0(-Po -m a + n a ) 9(p + m a - fx a ) | 1 

X \l + e /3(-po-m a +M 1 + e( 9(po+m K -Ma) + U[P ° + m ° ^ J 

-VaqXpoiPo - 2m a )e-^^- 2m ") {0(-p o ) 

f 6(m a - fig) fl(-m a + // a ) . 

X \ 1 + e /3(m a -^) - I + ^(-m.-^) ~ ~ W + °(\P\)- ( 60 ) 
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It is easily seen from eq. (pOj) that the additional term n^ =1 (0)Thcrmo be- 
comes zero at any temperature. 

Thus, from eqs. (57), (^) and fl59|) with the condition (^j), we find 
that eq. ( |56|) holds when the fermion massless limit is taken or when the 
chirally symmetric case (Nj~ = Nj = N) with massive (light) fermions is 
considered. In those cases, the following statement is allowed. According to 
the spontaneous magnetization, the magnetic part of the field strength has 
no any static mass and become long-range forces. If there is no spontaneous 
magnetization, the magnetic field becomes short range force due to the non- 
vanishing topological mass (k=k^0). 



7 Discussion 

In Sec. 6 we have shown that the spontaneous magnetization induces the 
fact that the static masses of magnetic fields vanish and the magnetic force 
becomes a long-range one. In spite of the direct breakdown of the Lorentz 
invariance by temperature effects, this relation between the magnetization 
and massless-ness of the magnetic fields is quite interesting, though the con- 
nection with the Nambu-Goldstone theorem has not been clarified yet. In 
addition, eq. (|56| ) is indeed consistent with the result from the thermody- 
namical potential approach. 

Now let us consider the reason why the spontaneous magnetization oc- 
curs even at high temperature. This result might look somewhat strange 
when we remember various symmetries which are spontaneously broken at 
low temperature and restored at high temperature; chiral symmetry, gauge 
symmetry for weak interaction, some physics of magnetization like in fer- 



romagnetic materials |Z3 and so on. At T = case, the viewpoint from 
observing the shift in zero-point energies of photons could explain the ap- 
pearance of the negative coefficient of the linear \B\ term to some extent ||. 
Similar approximate estimation at finite temperature can be also useful. In 
the finite temperature system, the free energies of photons should be con- 
sidered as well as the zero-point ener gies. If the re is no magnetic fields, the 
photons have the dynamic mass ~ J2Nuj 2 + k 2 [[L7[]. On the other hand, if 



the spontaneous magnetization occurs, the mass is shifted to ~ ^j2Nu 2 . The 
difference of the free energies of the photons between these cases is expressed 
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as 



AF boson = F boson^ m 2 = 2N J£ _ ^boson^ = 2N(J 2 + ^ 

= A£ z . p . + AF boson , (61) 
where AE z p ^ is deference of the photon zero-point energies and 



AF boson = 1 f OP e = (62) 

The rough estimation for AF hoson shows that there is the linear \B\ term 
with negative coefficients for small \B\ region, whereas the linear term in 
A-Ez.p. vanishes at high temperature limit. On the other hand, it can be 
seen that the positive contribution of B fields to Q f Q rmi amounts to the order 
0{\B\ 2 ) in the fermion massless limit. (When fermions are massive, there is 



the linear term with positive coefficient at T = 0. |10|, |2q]) These estimations 
are summarized in Appendix 1. Therefore we can realize to some extent that 
the spontaneous magnetization occurs even at finite temperature at least at 
small B cases. 

We have seen that the coefficient C{f3) becomes larger in the high temper- 
ature region (see fig. 3). We stress that this does not always mean that the 
spontaneous magnetic field, which is determined by observing the minimum 
of the potential, is an increasing function of T. Rather, it may be a decreas- 
ing function if the contribution of the higher orders of B to JT ad becomes 
dominant and it holds the value of the potential upper as temperature grows 
in the large B region. Then expectation value (B) would become smaller, 
even though C(f3) become larger as T grows. In this case, (B) goes to vanish 
effectively at high temperature limit. Such situation is very similar to the 



restoration of the parity violation analyzed by Klein-Kreisler et.al. [0. We, 
however, do not discuss any more the contribution of higher order of B to the 
thermodynamical potential Q(B,T) because our calculation here is basically 
due to the weak field approximation. 

When we consider this model in a realistic system such as a planar system 
in 3+1 dimensional world, the spontaneous magnetic field keeps to take non- 
zero values up to a certain critical temperature where a kind of the hopping 
to 3+1 dimension would occur and the system would be no longer regarded as 
a planar system. These characteristics may be interesting in the connection 
to the physics such as the high Tq superconductivity M. 
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8 Conclusion 



In this paper, we have investigated the spontaneous magnetization in 
QED 3 with a bare C-S term in the finite temperature system. At first, we 
have analyzed the thermodynamical potential Q as a function of temper- 
ature, chemical potential and magnetic field in the fermion one-loop level 
and the massless limit. Through these analyses, the consistency condition 
at finite temperature have been assumed, which gives the relation between 
/x, \B\ and j3. We have found that the negative coefficient of the linear B 
term is induced in Q by radiative correction. It means that the spontaneous 
magnetization can occur even at high temperature as well as at zero temper- 
ature. Secondly the spectrum of the physical photons has been investigated. 
At T = 0, the photons become massless by radiative effects so that they 
are regarded as the Nambu-Goldstone bosons for the spontaneous Lorentz 
symmetry breaking due to B ^ 0. We found that in our model, the topolog- 
ical mass vanish in the loop level even at finite temperature if the classical 
magnetic field exists. Since the magnetic force is a long-range one even at 
finite temperature unless there is a non-vanishing topological mass IIt p , we 
find that the following statement is allowed. According to the spontaneous 
magnetization, the magnetic part of the field strength has no any static mass 
and become long-range forces. This statement is consistent with the results 
by the thermodynamical potential approach above. Therefore we conclude 
that the spontaneous magnetic field, which is generated at T = 0, does not 
vanish even at finite temperature. This result is valid for the case of the 
fermion massless limit. 

As to the case with massive fermions, we have found that n topo vanishes 
in the chirally symmetric case (N + = N~ = N) putting q a = e > and 
m a = m ^ 0. However even in this case it remains unknown whether the 
spontaneous magnetization occurs at finite temperature in the model with 
massive fermions. We have not studied the details of the massive fermion 
cases by thermodynamical potential approach yet. This is a future prob- 
lem. These cases would be more important in connection with the fractional 
quantum Hall effects and high T^-superconductivity. 
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APPENDICES 



APPENDIX 1: 

The Hamiltonian operator, H = J2a H a , and the charge operator, Q 
Y^aQa, are given by 



^ = EE a lp a np^n + E E h lp h npU a n "EE W ». 



Qa 



"np^np^n 

n=0 P n=l p n=l p 

f Qa J2p{ a 0p a 0p ~\) \ 

\ Qa J2p(—bopbo p + §) J 

oo 

^Qa E E( a np a?1 P _ bnpbnp) , 
n=l P 



(A.1) 



> 0) 



• (A.2) 



The fermion sector of the thermodynamical potential (eq. fll^D) is ex- 
pressed in the classical level as 



Qfcrmi 
1J 



2 5 /% 2 



c (I EI^I 3/2 -E 



\q a B\ 



CXJ 

Eln(l + 



-0(uj n -q a fJ. a ) 



,n=l 



E In (l + e -^«+<^ a )) + In ( e /»^W2 + e - 



(A.3) 



n=l 



where we assumed the fermion massless limit (m a = 0). In the low temper- 
ature limit (/3 — > oo), the second term of r.h.s. vanish and the result in ref. 



is reproduced. In the case of (5 \q a B\ <C 1, Qq" 11 is calculated as |26 



Qfermi 
"0 



E 



+ 



Ql 



aS) + 0((3\q a B\ 
1 



dx 



Air [Jo 1 + ev^ 



+ o(\b\)\i4 + o0mI) 



(A.4) 



We see that there is no linear B term in f^Q 6 ™ 1 because of eq. (14). 

The estimation of AF is approximately performed as follows. Since 
the TliThermo rapidly becomes small for large |p| by the factor e~? la , we 
can employ the cut-off l~* e for the \p\ integral in AF boson , where l~^ e = 
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J2 IVall^al/ \VaQa u a\ ■ \B\ is defined in ref. ||. In the case of small l a 2 e / k 2 
(namely small \B\), we evaluate AF boson as 



1 rial i- e -PV^^I 
AF hoson = — - / dy In ' 



A-np Jo i_ e -Py/y+2Nwl+K* 

1 r l uve f 1 y + 2Nu 2 , ' 



47T/3 7o [2 y + 2A^ 2 + « 



p 

1 K 2 

"ttuei 



87T/3 2iVcj2 

B\ + 0{\B\ 2 ). (A.5) 



8tt In 2 1 

Thus we can see in this rough estimation that the linear \B\ term with a 
negative coefficient appears. 

APPENDIX 2: 

The fermion propagator in the magnetic field, S(x,y) = — i(0\Ti/j(x)ijj(y) |0), 
satisfies the relation 

S(x,y) qB< o = 
The Fourier transformation is given by 

S (P)= J -^^S Q {x-y), (A.7) 

where 

S(x, y) qB > = e-^y^+y^ 212 S (x - y) qB>0 . (A.8) 

The relation of the propagator between r] a = 1 fermions and r] a = — 1 ones is 

S (p) Va =-i = S (-po,p) Va=+ i. (A.9) 

The expression by the proper-time method is convenient in the calculation 
M. For the f] a = 1-fermions, the expression for v = is 



SoipWzt 1 = -H 2 / dr(cosr)- 1 e i ^- m T ~ tT T l tanr 
Jo 

x I (m cost + iposinr) / + ip l a\ — ip 2 er 2 ) (cost) -1 

+ (po cos r + im sin r) cr 3 } . (A. 10) 
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As to the case of v — 1, it takes the form 

Soipyzi 1 = Soipyzt 1 + f( P ), 

where f(p) is defined by f(p) = 27rie~^ 2/2 5(po — m)(I + as). 



(A.11) 



APPENDIX 3: 

Boson self-energies, IVs, in the weak field approxima- 
tion 

In our approximation, IVs are formally expressed as 
- 2 - d 3 k 



& I (2t^ 3 



m 



x 



n?=°(p,B) = 



-I + k (k - p ) + k 2 -p-k] ^ 

(k 2 — m 2 + it) \{k — p) 2 — m 2 + it] 
<->„• 2 „ j3 h. 



p 2 



d 3 k 
W) 



m. 



...aP 2 - i(p - 2k )(k 1 p 2 - k 2 pi) + q, B 2\ 
(k 2 — m 2 + it) [(k — p) 2 — m 2 + it] 
Aiql r d 3 k 



(A.12) 



(A.13) 



(p 2 ) 2 J (2k) 



x 



m; 



n™ =1 ( P ,B) = n a r°(p,B) 

^-nq 2 a r d 3 k m a + k 

J (2vfkl-k 2 -m 2 a + i> 



r 

; o(fco - Po) + k 2 -p-k^+p 2 [ml - k (k - p )} 
(k 2 — m 2 + it) \{k — p) 2 — m 2 + it] 

+0(B 2 ), (A.14) 



x 



{e-M 2 ^ 5(k -p -m) + e'^ 1 ' S(k + p - m)} ,(A.15) 

T-ry=0/ n\ 



nr =1 (p,5) = nr°(p,s) 

./ (2vr) 3 



p2 



X 



('Pi-^-(Pi + ^) fc i 5(A;o _ pQ _ ma)e -^ 2 ' 2 
k% — k 2 — m 2 a + ie 
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-g + (pi - tp 2 )fcl + (jgl + ^2^2 

(fc — p) 2 — m 2 a + ie 



6(k - m a )e 



-Pi* 



nr(p,B) = n-= u (p,B) 



(A.16) 
(A.17) 



The straightforward calculation of (|A.12[) ~ ( |A.17| ) reproduces the results in 
g |7j up to the order 0(B))] 



rr =0 



92 
8tt 



; ^-r^ I V^ + (1 - z ) Sl n 1 7F= 



y=0 



VaQa r ■ -1 

■— — v^sm 
4tt 



tt =0 
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where z = —4m 2 /p 2 and 



VTT^' 
^ (v^+(i-z) sin_1 ^ 



n; 



n 







1 



1 



«2 J_f 

2vr Po^ 1 P 2 + 2m a p P 2 - 2m a p 



TTai^=l 
ix 2 



nf =u - 



TTa!/=0 

11 2 



9a ?7a 



+ 



1 



27r Z 2 1 p 2 _|_ 2m a p p 2 - 2m a p 



(A.18) 
(A.19) 
(A.20) 



(A.21) 

(A.22) 
(A.23) 



APPENDIX 4: 

Iio, III and II2 in weak field approximation 

By using the replacement (pi]) and (p2|), the self-energies in the finite 
temperature system are formally obtained from eqs. (|A.12|) ~ ( |A.17|) . After 
performing the angular integral included in dk, the expression of HVs takes 
the following form. For the v — case, 



W=°(P3,P\B,(3) 



x 1 + 



q 2 1 ^ f°° fc 
"^2^5 ^ Jo fc| + A; 2 + m 2 

fc 2 + m 2 - p 2 - (k 3 - pa) (3fc 3 - Pa) 



/{(fc + |p1) 2 + m 2 + (fc 3 - p 3 ) 2 } {(fc - bl) 2 + m 2 + (fc 3 - p 3 ) 2 } 
+ 0{B 2 ), (A.24) 
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dx- 



+ 0(B 2 ), 



2tt (3 n ^ 00 Jo m 2 + (k 3 - p 3 x) 2 - x(l - x)^ 

(A.25) 



oo U 

dk- 



2/„2 ( A„2 i Q^2M.2 , (au ^ \ 1^1 i (az& „2^2 

Pi " 



p\k 2 - (4p 2 + 3^)fc| + (Ak 3 - p 3 )(p 2 + pl)p 3 - (Ap 2 - p 2 )m A 



sj{(k + \p\) 2 + m 2 + (k 3 - p 3 ) 2 } {(k - \p\) 2 + m 2 + (k 3 - p 3 ) 2 } 
+ C(5 2 ). (A.26) 



As to the case of v = 1, 



+ 



2irl 2 p 3 \p 2 + p 2 — i1p 3 m p 3 +fP + i1p 3 m 



m x C(e" 



-p 2 /(9a|B|)> 



(A.27) 



1 1 



+ 



2 

g 77 



+ 



2-7r/ 2 [ p 2 + p 2 — i2p 3 m p 2 + p 2 + i2p 3 m 

urHp3,f,B,p) = nr (P3,f,B,p). 



( e - P - 2 /MB|)) ; 



(A.28) 
(A.29) 



APPENDIX 5: Calculation of II 



1 Thermo 



(p = 0) 



The second term in r.h.s. in eq. (|58 ) is given by 



niThcrmoQo, B , (j , fx) = — I d\ 



o 2 



2iri fr ./ - x 



(27T) 



/ioc+fi—e I 
dk f(k ,X) _«( hn _, A , i ~ i„dk f(k ,\)\ , 



e -/3(fc -M) 4- 1 Jc 



(A.30) 
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where 



kg — k 2 — m 2 + ze 

~f + (gl - »P 2 )fcl + (jgl + P2)fc2 c i(fc - m) A c -Ptg 

(fco — Po) 2 ~ (k ~ P) 2 — m 2 + ie 
The straightforward calculation produces 

niThermo(P, B , (3 , fi) 



(A.31; 
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dcpe'^e 212 P0 (PO+2m a cos </>+p 2 sin 





Xe -P 2 l 2 e -P0{P0+2m a )l 2 



X 
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9(m a - /i a ) 



Po(Po - 2m a ) {0(-p o ) - 0(po - 2m a )} 

0(-m o + /ia) 



]_ _|_ e P(m a —IJ, a ) ]_ _|_ g/3(-m a +/^ a ) 
2tt 

o 



xe -F 2 ' 2 e -Po(PQ-2m a )« 2 



- 0(m a - /i a ) 



2i 2 -y/ po (po - 2m a ) (pi cos <^>+p2 sin </>) (A 32) 



To perform the integral for d(f) in eq. ( |A.32| ), we expand the integrants by 
\p\. Then the integrations become 



2tt 



)e =R</> e ±2Z 2 v /po(po±2m)(pi cos 0+p 2 f 



= ±2nl 2 ^p (p ± 2m) fa T ^2) + 0(p 2 ). (A.33) 
By using eq. ( |A.33| ), we finally obtain eq. QBOD from eq. ( |A.32| ). 



APPENDIX 6 

The formula to divide functions into the zero-temperature part and ad- 
ditional finite temperature part is W7 
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aE/W = o~ / dk f(k ) + —idk f(k ) 

1_ /-*°°+^ dk f(k ) _ J_ /-ioo+M-e dk f(k ) 

2ni J-ioo+^+e e^-ri + 1 2m J-ioc+n-e e-/ 3 ^"^^) + 1 ' 

(A.34) 

where /(&o) is a function without having singularities on imaginary axis and 
contour C is defined in Fig. 4. 

Euler-Mclaurin's formula is given by 

TO | ra+nh | 

^ g(a + hk) = - g(x)dx + - {g(a) + g(a + nh)} 

k=o nJa 1 

(A.35) 

where 5 r = 2(2r)!/(2 2r - l)7r- 2r En=i( 2n + f)~ 2r (the Bernoulli numbers). 
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c 


-/ oo 
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Figure 4: The contour C for the integral in eq. (A.34). 
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